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Abstract Perturbations of satellite orbits in the gravitational field of a body 
with a mass monopole and arbitrary spin multipole moments are considered for 
an axisymmetric and stationary situation. Periodic and secular effects caused by 
the central gravitomagnetic field are derived by a first order perturbation theory. 
For a central spin-dipole field these results reduce to the well known Lense-Thirring 
effects. 
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1 Introduction 

As is well known Einstein’s theory of gravity leads to so-called gravitomagnetic 
effects, i.e., gravitational effects caused by moving or rotating masses. Such effects, 
e.g., related with a rotating Earth, have been studied intensively in the past. To 
be mentioned here is the Lense-Thirring effect that describes perturbations of 
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a satellite orbit in the gravitomagnetic dipole fie ld of a rotating body. It was 


first described by Lense and Ttiirring (1918, see iMas hhoon et al. (Il984ll for a 
translation) after extensive groundwork by Einstein. ICiufolini and Pavlid ( 2004ll 
found the effect to be confirmed within a 10 % accnracy by a detailed analysis of 
the orbits of the two L AGEOS satellites. The resnlts, however, are still snbject to 
an ongoing debate (see llorio et ^ (l201lh ~l. 

Closely related to the Lense-Thirring effect is the precession a torque-free gyro¬ 
scope experiences due to gravitomagnetism. T his Pugh-Sch iff e ffect was pro posed 
to be an alternative test for frame- dragging (see Pugh ( 1959ll an d lSchifj(ll96nli i. Al- 
though unexpected problems arose l^eritt et ah (|2011 1 measured the effect within 
an accuracy of 13 %. 

The above mentioned effects solely take the spin dipole into account. Natu¬ 
rally one is also interested in the influence of higher spin multipole moments. Such 
higher multipole moments come into play via multipole expansions of the so-called 
gravitoelectric potential 0 (a generalisation of the Newtonian potential U which 
is often also denoted by w) and the gravitomagnetic vector potential w. Both are 
used in order to parametrise the metric tensor in the first post-Newtonian approx¬ 
imation. For more details on the origin an d form of th e metr ic ten sor in the first 
post-N ewtonian approximation please see ISoffel et al.l ( 2003h and iDamour et ^ 
(Il99lll . 

In this paper w will be our primary concern. As for the gravitoelectric potential 
(j) we will just use a post-Newtonian mass monopole, i.e., (j) = GM/r, where G is the 
Gravitational constant, M the Blanchet-Damour mass of the central body and r = 
||x|| the distance to its centre of mass. The gravitomagnetic potential w is induced 
by a matter current density cr (cr^ := j where is the body’s energy- 
momentum tensor and c the speed of light). In the stationary case and outside a 
coordinate sphere S that fully covers the energy-momentu m tensor of the central 
body w admits a multipole expansion of the form (e.g. iBlanchet and Damouil 
(Il989|) ) 


w = Wfeefc 


= -gE 




^kabJbL — l^aL—l 


®/c 


with the spin multipole moments 


/ 

K3 


• / ^abKki^L — l^a^b^ 


( 1 ) 


( 2 ) 


In these equations Cabc denotes the fully antisymmetric three dimensional Levi- 
Civita symbol for which we shall use £123 := +1. The vectors (e^), k = 1,2,3 
stand for the canonical basis of the i.e., ei = (1,0,0)^ etc. The spin-moments 
Jl are Cartesian STF (Symmetric and Trace-free) tensors, where L is a Cartesian 
multi-index, L = (A:i ,... ,ki) of Z Cartesian indices, each taking the values 1, 2, 3. 
Symmetric refers to the symmetry with respect to all I indices while trace-free 
means that every contraction between two arbitray indic es va nishes. For more 
i nform ation on STF-tensors please also see iThornel ( 198011 and iHartmann et al.l 

(I1994I) . 


A summation over two equal (dummy) indices is implied automatically (e.g., 
AaBa := AiBi + A 2 B 2 + A 3 B 3 ) and dab '■= /{dx°‘dx^) etc. The hat on top of a 
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symbol indicates the STF-part. Angle brackets indicate the STF-part with respect 
to the indices enclosed. 

Note that for the hrst post-Newtonian metric the spin multipole moments have 
to be defined to Newtonian order only. The spin-dipole moment, J, agrees with 
the usual intrinsic angular momentum of the body. 

In this paper we will study a central axisymmetric body rotating uniformly 
about its symmetry axis. Perturbations of satellite orbits caused by the correspond¬ 
ing gravitomagnetic field induced by spin multipole moments of arbitrary order 
will be considered by me ans of perturbation theo ry. Such spin multipole mo ments 
have been considered bv iTevssandieij ( 1978jl and IPanhans and Soffell (|2014|). Ear¬ 
l ier ap plications of the mentioned expansions ca n be found inizimbres and Letelieil 
(1200811 for the precession of a gyroscope and in llori 3 (l200lh where an alternative 
derivation of the Lense-Thirring effect is given. 

Here we shall also use a coordinate system with an axial symmetry with re¬ 
spect to the z-axis, i.e., the spin vector of the body points into z-direction. This 
simplihes the form of the used multipole expansion. On the other hand this might 
be a shortcoming because of the loss of generality due to unknown transforma¬ 
tion behaviour of the spin moments. A treatment of a general spin orientation 
can be found in llori 3 (I 2012 D where a general spin-spin interaction was taken into 
account as well. A discussio n about mix ed effects connected with an arbitrary spin 
orientation can be found in llori 3 (|2015[) . 


2 Equations of motion and STW-decomposition 


Post-Newtonian satellite equat i ons of motion have been studied in detail in the 
literature (e.g., IPamour et al.l ( 1994h l. Considering only a single isolated cen- 
tral body the coord inate acceleration of a satellite as (e.g., equation (3.4) of 
IPamour et al.l ( 1994ll l has the form 


as = V(j} + ^ 1^—4(?!)V(?!) — v(30 -|- 4(vV)(?!)) -|- 4w -|- -|- v x B | (3) 


where the gravitomagnetic field B is dehned by 


B := -4V X w . (4) 

Since our central interest are the perturbations induced by some stationary grav¬ 
itomagnetic field, i.e., w = 0, we will simplify this equation to 


as — V0 + aper 


( 5 ) 


with the perturbing acceleration 

aper := 4jV X B . 


(6) 


Here, v is the satellite’s coordinate velocity. Effects connected with a varia t ion in 
tim e of the spin vector and orders of magnitude were studied in iBini et al.l ( 2008ll 
and llori3 ( 2002ll . For satellite orbits around the Earth such time dependencies can 


be neglected. 
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As mentioned in the beginning we will take 4> = GM/r so the gravitoelectric 
part of the potential leads to unperturbed Keplerian orbits. The satellite orbit 
will be described by the usual set of orbital elements (a, e,/, 17, cu, Mq) and the 
perturbations by means of a Gaussian STW-perturbation theory of first order 
with S', T, and W being the scalar products of aper with n for S, k for W and 
n X k for T. These vectors are defined by n := x/r, r := ||a;|| and k := C/G, 
C:=xxv, G:= ||C||. 

An explicit calculation yields 


S = 


B C 


rc^ 




C 


Gc2 


( 7 ) 


In a next step we want to calculate B explicitly. Using o a calculation of the 
B-field reveals 


B = 4G^ 




JlVOl 


J- 


( 8 ) 


The use of spherical coordinates has proven to be of advantage and so we will use 
spherical spin multipole moments 


^im ■■= J d^x p(x X 

as introduced bv IPanhans and Soffell ( 2014ll . Here 




stands for spherical harmonics and 


.— I 21 + 1 (I - m)\ imtf, 


All {I + m)! 


Plm{\) ■■ = 


i-iy 


■{l-A )2 




PlmW 




2H\ ^ ^ 

are associated Legendre functions. Legendre polynomials are denoted as 

Pi{\) := Pio{X). 


( 9 ) 


( 10 ) 


( 11 ) 


( 12 ) 


In the following the arguments of the these functions will be the polar angle 
while we plug in cos{0) for A where 0 is the azimuth angle. With this convention 
we will suppress arguments when using these functions. 

The spherical spin multipoles are connected with via 


Jl = 


Airil — 1 )! 

(2Z + 1)!! ^ , 

m= — L 




Written in terms of spherical harmonics B becomes 

CXD I 




47r 


^lm=-l (2^ + l)(^ + l) 


^ / 4 


VYlm -{1 + l)Ylr 




} 


(13) 


(14) 
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The assumption of an axial symmetry implies ^im = Under this assumption 

expression m reduces to 


B = -4G^ 
^ = 1 


V (2^ + 1)(^ + 1) 




sin(g)P/ 


o/-|-2 


eg + {I + l)Pi- 


. 1+3 


(15) 


with 

eg = cos(0) cos{(j))ei + cos(0) sin((())e 2 — sin( 0 )e 3 (16) 

while Pi is the first derivative of Pi with respect to its variable. 

As consequence of this form of B we get for S and W (T is given by S through 

©) 


S = 

w = 


8^CG 


E 


& 


c2 ^ y2r+T(^ + i)H+3 


cos{I)Pl, 


8y/^CG 


E 




c2 ^ y2r+T(^ + l)r'+3 


{I + l)Pl + 


re sin(t/) cos(u) sin(7) 
P 



(17) 


(18) 


In this equation v notes the true anomaly and we used u := v+ui and p := a(l — e^). 
In the next section the perturbations of the orbital elements will be discussed. 


3 Discussion of the single orbital elements 

3.1 Semi-major axis 

The differential equation for the semi-major axis reads 

Sesm(iy)+T^) . 
We make use of the second equation in © and get 

2 7 . , ^ X- VpNj 

a = — , esiniiy) -- b. 

n^/D^ V ^ ^ G rJ 

In the last step we furthermore apply 

tG . , . 

X • V = —esm(i^}, 

P 

an equation we will also need for the eccentricity, and find 

a = 0. 



(19) 


( 20 ) 


( 21 ) 


( 22 ) 


So for this special form of the perturbation acceleration the semi-major axis will 
not change over time. 
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3.2 Eccentricity 

We are using again the second relation of © and the expression for x ■ v together 
with equations known from the description with orbital elements, 


T P 

cos(E) = -(e-\-cos(iy)), - = 1 + ecos(i^), 

p r 


(23) 


and find 

. _ 8y/nGcos{I) 


E- 




j sin(i/)f’/(sin(u)sin(/)). (24) 


c2 ^/2^hT(^ + l)a*+2 

The last three steps are expressing the Pi with via 

i-i 

Pi = E = 


k=0 


2fc + 1 k even, I odd or k odd, I even 


0 


otherwise 


(25) 


(a consequence of no. 8.915/2. in iGradshtevn and RvzhikI ( 200l1l') followed b y 
rewriting the Pk in terms of complex inclination functions i/taf, fsee lKau i3 (|l96lll i 


f’fe(sin(M)sin(/)) = ^ Efc 06 (/)e*^'' 


(26) 


6=0 


and finishing the series of conversions by eliminati ng the implicit time dependency 
throngh i/ by using Hansen coefficients fe.g.. lBranham Jrl ( 1990lH so we end 

with 


8y/n cos{I)G 


E 


l — l k 




EE 


1=1 

4 E (.v,-‘-»-“«-.v. 


-l-l,k-2b- 


l\ ii 


M 


(27) 


s= —oo 


In a first order perturbation theory we just have to integrate this expression with 
respect to t which is fairly easy seeing the trivial time dependency M{t) = nt+Mo- 
One has to be careful with the term for s = 0 though. The result is 

l-l k 


Ae = 


8^/n cos{I)G 


E 


6 


r 1 




— 


EE NikFkOb{I)P'‘^^’‘ 


+ 


s= — oc> 
s^O 

i 


-l-l,k-2b-\-l v-l-l,k-2b- 

^0 




(28) 


The secular perturbations of e vanish because of 


Aesec{t) = 


8tG^/^T cos {I)^i 


E 


i 

= 0 


^ v/2m(/ + l)a*+2c2 

fe=even 


2i 


\r — 1 — /,1 T/’-1 — I, — 1 

^0 “ ^0 


( 29 ) 


i_ J -^rn.m ■\rn. — m 

where we used Aq = Aq 
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3.3 Inclination 

The discussion of the inclination I will follow the same pattern as before but gets 
slightly more complicated because of the appearance of W rather than S. The 
starting point is the perturbation equation 

(I + \)Pi + -esm(iy) sin(/) cos(u)P[ . 
P 

(30) 

This time we pass on expanding P( but use 

cos(m) sin(7)T’/(sin(M) sin(7)) = ^^(sin(M) sin(7)) (31) 

instead. We apply the other conversions as before and hnd 
f ^ Ay/nG ^ 

+ 1(^ + l)a'+ip Vn/ 

i 

Fiob (^e*“ + e“*“^ _|_ i)(i + ecos(t^)) + ie(l - 2b) sin(i/)j 

6=0 

(32) 


I = 


Sy/^G 


E 


Cl 


V2l + 1(1 + l)r'+2 


cos(m) 




E 


c2 ^ 72 r+T(Z + l)a'+ip 


0 ""E^-e 


iuj{l-2b) 


b=0 


I [(^ + l)e' 


+ ^{21 - 2b+l)e 
+ ^{2l-2b+l)e 


+ ^{2b + l)F 




^ I 


Y _^ p ru(l-2b) 

^y/2l + T(Z + l)a'+2(l-e2)c2^^ 

OO ( 

^ \ [(/+ +-(2^ - 26+ 

S=-CX) I ^ ^ 

+ |(26+ l)X7^“'’*“^'’j + [(/ + 

+ |(2/ - 26+ 1)X7^“*’'“^'’ + |(26+ l)X7^“*’'“^'’“^je“*‘^|e'^ 


(33) 
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An integration yields the fairly long expression 


= E 


4G0rf; 


I- VWTTil + l)a'+2(l - e2)c2 ^ 


E ^«0b' 


^iui{l-2b) 


l=l 


b=0 


oo /• 

?J[ 


(i + i)x; 


+ -(2/-26+l)X; 


s= —oo 
s^O 




— 1 — 1,1 —2b icj 


e*"" + (/ + 1)X, 




+ |( 2 /- 26 +l)X 7 ^“'’*“^^ + |( 26 +l)X 7 ^“'’*“^'’“^' 


le-"|5^ 

J J tns 


{[ 


+ < \{l + l)Xo 


+ -(2/ - 26 + 1 )Xq 


-l-(.i-2b+2 


+ I (26+ e*"" + [(/ + 1)X, 


-l-l,l-2b-l 

0 


+ -(2l-2b+l)X- 


' + -{2b+l)X- 




(34) 


As for the secular perturbation we find 


A/sec (f) — 


E 


c2 ^ V^TTTil + _ g2) 

^=odd 

[(/ + + 2)X-i-*--i] 

+Fioi^ [(/ + + |(/ + 2)X-^-^’^ + } 


E- 


fi 


4y^Gi 

c2 ^ 72 r+T(i + l)a'+2(l-e2) 

(/ + l) + eXo-i-'-i) 


= 0 , 


(35) 


which vanishes again due to F^^i+i + F^^i^ = 0 for odd numbers 1. 


3.4 Argument of the ascending node 

Luckily the above argumentation can be adopted for the calculation of Af2 with 
some minor changes since the differential equations for f2 and I just differ by 
the appearance of sin(u) rather than cos(m) and an additional sin(/) in the de¬ 
nominator which has no influence on the calculation at all. So the result for AO 
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reads 


Anit) = J2-. 




sin(/)^/2^^■(/ + l)a'+2(l - e2)c2 ^ 


Fwbe 


u{l-2b) 


-l-l,l-2b+2 


^ OO f 

S= —CXD ^ 

s#0 

+ |(26 + - [{/ + 


+ |(2Z- 26+ -L - 


+ 2 


-l-Ll-2b-2'] -ioj 1 e 


AsM 

n 

J J ins 


+ I [(Z + + |(2Z - 26 + l)Xo 


\-l,l-2b-\-2 


+ -(26+l)Xo- 


— \ — l,l — 2b~\ i<jj 


A - (Z + 1)X, 




+ |(2Z - 26 + V)X-^-^'^-'^^ + + 1 )^ 0 “'“'’'“^^“^^ ■“ 




(36) 

This has noticeable consequences in particular for the secular perturbations which 
read 


2\f2sec(Z) — 


8y/nGt 


E 


sin(/) V2Z + la'+2 

l=odd 


6 


{^10^^) 


X, 


- 2 - 1,0 

0 


S-y/TrCrt V 

c2sin(/) ^ y/2r+Ta'+2 (1 - e2)'+2 




t(f) 


n.=0 


/=odd 

e\2n ^2n\ / Z 
n / i 2n 


where we used 


a 2\\^n,m -v^n+ljm . ^ ( ^n+l,m+l , -v'n+l,m— 

— e jAg — Ag + - (^Ag + Ag ^ 


(37) 


(38) 


with s = m = 0,n = —l — 2 and 




_ I _ y (A' 

(l-ey-i to 


e\2b^26 + m\/ n—1 
26 + m 


(39) 


for n G N, m G Z. 

If onl y the spin-dipole i s kep t we recover the well known result of Lense and 
Thirring ( Mashhoon et all ( 198^ 1: 


. pi=i 8v^Gt __ c^tjr f 

- c2sin(/)^„3(i_e2)f'^^ c2a3n(l-e2)fV 3 


^int. (40) 
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Inserting the Lense-Thirring expression is obtained. 

Because of its relevance we want to apply this formula to the LAGEOS 2 
satellite. We will treat Earth as a homogeneous body. Since the quadrupole term 
does not contribute to the secular perturbations at all the dipole needs to be 
compared to the octupole term. Evaluating formula dsZl) gives 


Ai=3 _ o + je^) 

V 7c2a5(l-e2)i 


^4 cos^ 




+ 4 cos 




(41) 


for the secular drift rate caused by the spin octupole. For the model of a homoge¬ 
neous Earth one finds 

(42) 


Qi-p = 0.02—. 
yr 


Since our model of a homogeneous Earth overestimates the value of the spin oc¬ 
tupole moment the actual value for the drift rate is even smalle r than the calculated 
numb er. If one compares this value to the spin dipole term ( Ciufolini and Pavlid 

(l200l ). 

(43) 


n' = l _ qi c;"13,S 

^ ^sec — 31.5 , 

yr 


it becomes obvious that the contribution of the spin octupole and all higher spin 
multipole moments can be neglected for Earth’s satellite orbits at present. 


3.5 Argument of periapsis 


Most of the work for this orbital element is done by now because the result from 
[Q can be used. However, in the differential equation for ui appears an additional 
term Wadd which needs to be studied. So we calculate 


^add • — 




— S cos(i^) + T sin(r') 


1 + 


P\ 


C_ 
fie 

8y/npG cos(/) 


f 2 / ^ \ 2 

cosM -\ —esin (t^) + - esin (i/) 

P \P 


E 


ec2 ^ ^/2r+T(/ -h l)H+3 


/ \ 2 \ *-l k 

cos(i^) +-esin^(i^)-h f - j esin^(z^) EE NikFkob{I)e'^^^ 
^ '^7'/ J 


(44) 
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8y/npG cos(/) 


E 


ec2 + 


EE 


OO f ^ 

l( 


s= — oc> 


f v-l-3,k-2b-\-l I v-l-3,k-2b- 
1 « “r s 


0 


+ 


2(1 - e2) 

e 

2(l-e2)2 


-l-2,k-2b 




-5 IX 


-/-2,fc-26+2 , v-l-2,k-2b-2 

] -^s 


) 

y-l-l,k-2b _ ^ Y-l-l,k-2b+2 , y-l-l,k-2b-2'\ 


(45) 


We integrate (1451) and find 

= -y ^ff^cos(/)e L_ y y ^ p, ^jyuik- 2 b) 

)^ec2^/2ZTT(/ + l)a'+3 


OO ^ ^ 

i:{l( 


—I — 3,k — 26-1-1 — I — 3,k — 26 — 


‘) 


s^O 


2(1-e2) 
e 


j^-l-2,k-2b _ ^ ^^-l-2,k-2b+2 j^-i-2,fc-26-2 ' 


2 (l-e2)2 

e 




l(^ 


-Z-l,fc-26-H2 I -v^-Z-l,fc-26-2 

“r s 


/ J ms 


■I/' — l — 3,fc — 26-|-l I — I — 3,k — 26—1 

"r 


+ 


+ 


2(1 - e2) E° 


-l-2,k-2b 


X, 


2(l-e2)2 


( 

K- 


-/-2,fc-26-H2 , v-l-2,k-2b-2 

+ A, 


0 


0 


- ^ —l,fc —26 / T2- L — l.,K, — ZO-\-Z I -^- 

^0 “ o + ^0 


^-l-l,k-2b-\-2 . ^-l-l,k-2b-2\ 


)] }| 


(46) 


So using the result from 13.41 we find for the overall perturbation 
Aoj{t) = — COs(/)Z\f7(f) + Z\(2;add(fl) 

with Af2{t) from (1361) and Zitdadd(^) from (1461) . The secular perturbations are given 
by 


ec2^/2ITT(^ + l)a*+3 E 


ZiWsec (f) = - cos /Z\12sec (f) - ^ C0s(7)fCi 

e{ 

3.1 , V 


V—^I 

An -h ■ 


f -y- I — 2,0 -y — I — 2 

' ^0 “ ^0 


2(1-e2) 




=0 

fc=even 

■y—^ —1^0 y — I — 

^0 “ ^0 


■") 


2(l-e2)2 


( 47 ) 
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Again, for the spin-dipole the result agrees with the classical expression of 
Lense and Thirring: 


1=1 

^sec 


- cos(/)Z\J7Lc^ - 

- 3cos(/)Z\J7'=^ 


8^/^^Gp cos{I) e ^ 

ec2 ^a4(l-e2)f 


(48) 


4 Summary 

The aim of the paper was to investigate the influence of a central gravity field with 
a mass monopole and arbitrary spin multipole moments on satellite orbits given 
a stationary axisymmetric setting. In order to simplify the form of the multipole 
moments a coordinate system was chosen for which the spin vector pointed into 
the z-direction. We found that a = 0 holds in general (and not just in first order 
perturbation theory) and that e and I experience perturbations which have no 
secular contributions. For odd numbers I there are secular perturbations for 
and CO which yield for the spin-dipole, ^ = 1, the well known results by Lense 
and Thirring. In the case of the LAGEOS 2 satellite we calculated the additional 
secular drift due to Earth’s spin octupole and found a negligible small number. 
Eurther Physical implications and orders of magnitude will be discussed elsewhere. 
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